The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

超临界压力CO2在水平圆管内流动传热数值分析

采用SST k-w低雷诺数湍流模型对加热条件下超临界压力CO2在内径di=22.14 mm,加热长度Lh=2440 mm水平圆管内三维稳态流动与传热特性进行了数值计算.通过超临界CO2在水平圆管内的流动传热实验数据验证了数值模型的可靠性和准确性.首先,研究了超临界压力CO2在水平圆管内的流动传热特点,基于超临界CO2在类临界温度Tpc处发生类液-类气“相变”的假设,揭示了水平圆管顶母线和底母线区域不同的流动传热行为.然后,分析了热流密度qw和质量流速G对水平圆管内超临界压力CO2流动换热的影响,通过获取流体域内的物性分布、速度分布和湍流分布等详细信息,重点解释了不同热流密度qw和质量流速G下顶母线内壁温度Tw,i分布产生差异的传热机理,分析结果确定了类气膜厚度d、类气膜性质、轴向速度u和湍动能k是影响顶母线壁温分布差异的主要因素.研究结果可以为超临界压力CO2换热装置的优化设计和安全运行提供理论指导.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

A modern approach to the solution and analysis of the Simha-Somcynsky model for the description of pressure-volume-temperature properties of polymer fluids

We revisit the Simha-Somcynsky model of polymer fluids with the purpose of developing novel theoretical and computational approaches to simplify and speed up its solution as well as the fitting of experimental data, and decrease its level of mathematical complexity. We report a novel method that allows us to solve one of the two equations of the model exactly, thus putting the level of mathematical difficulty on a par with the one of other models for polymer fluids. Moreover, we describe a computational algorithm capable of fitting all five parameters of the model in an unbiased way. The results obtained reproduce literature results and fit experimental pressure-volume-temperature and solubility parameter data for three polymers very accurately. Moreover, the new techniques allow for the investigation of the model at very low temperatures. Unexpectedly, the model predicts behaviors that could be interpreted as a glass transition, as routinely observed in dilatometry and differential scanning calorimetry, and a glass phase. We compared the predicted and experimental T g’s for cis poly(1,4-butadiene) and found an excellent quantitative agreement.     

浮升力和流动加速对超临界CO2管内流动传热影响

基于单相流体的概念,超临界流体的异常传热行为已经被研究很多年了,但是关于其流动传热机理仍没有统一的认识.本文通过理论分析和实验研究了超临界二氧化碳在竖直管内向上流动过程中,浮升力和流动加速效应对其流动结构和传热过程的影响.结果表明,没有确凿的实验证据表明超临界流体的异常传热行为是浮升力和流动加速直接导致的,存在的估计浮升力和流动加速效应准则均是在常物性流体的基础上,做了大量假设得出的,不同的研究者采用浮升力和流动加速准则分析超临界流体的传热恶化得出的结论不一致.最后,基于拟沸腾理论分析超临界流体的传热恶化过程,提出超临界沸腾数区分了超临界流体正常传热与恶化传热的转换边界,为超临界流体流动传热研究提供新思路,超临界沸腾数对建立用于不同技术的超临界流体动力循环的最佳运行条件具有重要意义.     

超临界水驱超稠油提高采收率热物理特性研究

超稠油资源的高效开发对于提高石油供给和保障我国石油安全意义重大。针对深层超稠油资源由于原始地层压力高、原油黏度大导致常规蒸汽驱无法有效开发的问题,本文提出超临界水驱开发深层超稠油的新思路。本文首先研制了超临界水驱油提高采收率实验平台,研发的管式填砂岩心模型能够模拟岩心升温和驱替的同步过程,然后开展了超临界水驱、蒸汽驱和热水驱对比实验研究,实验结果表明,相对于蒸汽驱,超临界水驱能显著提高采收率并具有更高的热效率,25MPa、400℃超临界水驱鲁克沁超稠油的采收率达到97.07%;获得了超临界水驱过程的温度场和驱替压差变化规律,发现了超临界水超覆现象,与蒸汽超覆相比,超临界水超覆发生晚且持续短,可扩大波及范围,提高采收率。     

铅基堆超临界CO2复合循环发电系统热力学分析

本文对比了再压缩超临界CO2 (S-CO2)循环、蒸汽朗肯循环、He布雷顿循环分别应用于铅基堆的最优热学性能,明确了S-CO2循环与铅基堆结合较传统循环的热力学优势。为进一步提高再压缩S-CO2循环的效率,以跨临界CO2 (T-CO2)循环为底循环构建了再压缩S-CO2/T-CO2复合循环,探讨了不同顶循环透平入口温度、压力和压缩机入口温度条件下系统性能的变化规律,对比了S-CO2/T-CO2复合循环和S-CO2循环的热学性能。结果表明:铅基堆再压缩S-CO2循环发电系统较传统循环形式具有更高的热效率;构建的S-CO2/T-CO2复合循环能够有效提高S-CO2循环的效率,在所研究参数范围内,S-CO2/T-CO2复合循环的热效率和效率比S-CO2循环分别最大可提高约4.8%和8.3%;再压缩S-CO2循环和S-CO2/T-CO2复合循环热学性能随顶循环关键参数变化规律具有一致性。     

Biocatalytic Feedback‐Controlled Non‐Newtonian Fluids

Non‐Newtonian fluids are ubiquitous in daily life and industrial applications. Herein, we report an intelligent fluidic system integrating two distinct non‐Newtonian rheological properties mediated by an autocatalytic enzyme reaction. Associative polyelectrolytes bearing a small amount of ionic and alkyl groups are engineered: by carefully balancing the charge density and the hydrophobic effect, the polymer solutions demonstrate a unique shear thickening property at low pH while shear thinning at high pH. The urea‐urease clock reaction is utilized to program a feedback‐induced pH change, leading to a strong upturn of the nonlinear viscoelastic properties. As long as the chemical fuel is supplied, two distinct non‐Newtonian states can be achieved with a tunable lifetime span. As a proof of concept, we demonstrate how the physical energy‐driven nonequilibrium properties can be manipulated by a chemical‐fueled process.     

超临界水中硫酸钠结晶动力学的分子动力学模拟

在超临界水反应器中，硫酸钠是易造成堵塞的一种常见无机盐，研究其结晶动力学对于防盐沉积反应器的设计具有重要意义.本文采用LAMMPS分子动力学模拟软件研究硫酸钠在超临界水中的微观结晶过程，其中水分子采用SPC/E模型，离子-离子、离子-水分子相互作用采用Coulumb和Lennard-Jones联合势能函数.结果表明：水对离子的静电屏蔽作用随温度升高而增强、随密度减小而减弱；增大超临界水的温度和密度有利于离子扩散，进而促进离子相互碰撞、成核；在模拟的超临界水参数范围内，其成核速率的数量级为10²⁹cm^-3·s^-1，生长速率为(19.8~25.8) m·s^-1.